Observational constraints from SNe la and Gamma-Ray Bursts on a clumpy universe 
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The luminosity distance describing the effect of local inhomogeneities in the propagation of light 
proposed by Zeldovich-Kantowski-Dyer-Roeder (ZKDR) is tested with two probes for two distinct 
ranges of redshifts: supernovae la (SNe la) in 0.015 < z < 1.414 and gamma-ray bursts (GRBs) 
in 1.547 < z < 3.57. Our analysis is performed by a Markov Chain Monte Carlo (MCMC) code 
that allows us to constrain the matter density parameter f2 m as well as the smoothness parameter a 
that measures the inhomogeneous-homogeneous rate of the cosmic fluid in a fiat ACDM model. The 
obtained best fits are (fi m = 0.285l£ u ^, a = 0.856^;^) from SNe la and (Q m = 0. 259±° ° ll, a = 
0.587^0 202) from GRBs, while from the joint analysis the best fits are (f2 m = 0.284to'o20> a — 
0.685^o'i7i) with a x? c d = 0.975. The value of the smoothness parameter a indicates a clumped 
universe however it does not have an impact on the amount of dark energy (cosmological constant) 
needed to fit observations. This result may be an indication that the Dyer-Roeder approximation 
does not describe in a precise form the effects of dumpiness in the expansion of the universe. 
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I. INTRODUCTION 

One of the most important problems in modern cos- 
mology is the coincidence in the order of magnitude be- 
tween the inferred energy density of the source driving 
the acceleration and the energy density of matter. Since 
the most evident qualitative change in the late universe 
is the formation of nonlinear structures, one cannot help 
noticing that cosmic acceleration and the formation of 
large structure occur, coincidently, in near epochs. If the 
recent formation of large structures could be identified as 
the driver for cosmic acceleration, the coincidence would 
be explained in a natural way. In fact there are theoreti- 
cal results that prove that gravitational collapse can lead 
acceleration, see for example pQ. Moreover, in Refs. [5] 
and [3] , it is argued that the influence of inhomogeneities 
on the effective evolution history of the universe is en- 
coded in backreaction terms, and that these terms can 
mimic dark energy; in this context see also Ref. [3] for a 
review on different approaches to the dynamical aspects 
of backreaction effects in recent eras. 

It has also been addressed by several authors (see for 
instance Refs. [3HS]) that the process of averaging and 
solving the Einstein equations is non-commutative and 
that maybe if we do it carefully, the accelerated expan- 
sion lately observed can be explained without the need 
of dark energy. 

On the other hand, most cosmological observations 
probe quantities related to light propagation such as 
the redshift, the angular diameter distance, the luminos- 
ity distance or image distortion. The prevailing model 
ACDM assumes that this propagation occurs through a 
homogeneous perfect fluid; but on scales smaller than 100 
/i^ 1 Mpc it is not reliable to consider that the universe can 
be modeled by a FRW cosmology. 

The incorporation to the FRW model of local inhomo- 
geneities have been addressed by several authors since 
the sixties, namely, works by Zel'dovich [7j, Bertotti [5], 



Kantowski .9] and, more recently, by Dyer and Roeder 
[TO] , again by Kantowski [IT] and by Demianski [12] . The 
general idea is to consider the propagation of light rays in 
a fluid that may clump. The proportion of clumped mat- 
ter respect to the homogeneous fluid is measured by the 
smoothness parameter. Departing from the Sachs expres- 
sion for a congruence of null geodesies, neglecting twist 
and shear effects, and including in the matter content the 
smoothness parameter, they arrive to an equation for the 
angular diameter distance, that can be connected, via the 
Etherington relation, to the observable luminosity dis- 
tance, resulting then the so called Zeldovich-Kantowski- 
Dyer-Roeder (ZKDR) luminosity distance. 

In this paper we consider observations of SNe la in 
the range of 0.015 < z < 1.414 [13] and GRBs in the 
range of 1.547 < z < 3.57 Q3] to test the ZKDR distance 
and constrain the smoothness or dumpiness parameter. 
For our analysis we employ the Markov Chain Monte 
Carlo (MCMC) technique, it leads us to obtain better 
constraints than previous analysis [15] . 

Furthermore, the GRBs sample is a very improved one 
with a significative diminished dispersion in the data 
[14] , From the joint analysis, SNe la and GRBs, the 
best fits are (0 TO = 0.284±°;°^,a = 0.685^;^) with a 
Xrcd = 0.975, indicating a clumped universe. Despite this 
fact, the amount of the dark energy (cosmological con- 
stant) needed to fit observations is not less than in the 
homogeneous FRW model. This may be an indication 
that the Dyer-Roeder approach is missing something in 
its physical description and should be revised, see Ref. 
[16) and references therein for a further discussion. 

In the literature the ZKDR luminosity distance has 
been probed also with Hubble parameter measurements. 
Similarly, we include this test using the data reported 
in Ref. [17], obtaining from the joint analysis (SNe la, 
GRBs and Hubble) the best fits: n m = 0.275^;^ and 

a = 0.82lt°;^° with a xLi = °- 974 that improves the 
previous analysis reported in Ref. |18) . 
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The paper is organized as follows: in Sec. II we briefly 
sketch how to obtain the ZKDR luminosity distance from 
the null geodesic congruence in a clumped universe, using 
an exact solution of a hypergeometric equation; in Sec. 
Ill the observational data and statistical method are de- 
scribed. In Sec. IV we discuss the results and finally in 
the last Section concluding remarks are given. 



II. THE ZKDR LUMINOSITY-DISTANCE 

Let us consider the light propagation in the geometric 
optics approximation. The two optical scalars, expan- 
sion and shear, 9 and a, respectively, satisfy the Sachs 
propagation equations [15] , 
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2 = - l -R a pk a k? , (1) 
Ic^sff^k^k 6 , (2) 



where a dot denotes differentiation with respect to the 
affine parameter A and fc M is the null vector field tan- 
gent to the light ray. Notice that since we shall address 
the RW geometry, the Weyl tensor is zero and the shear 
vanishes. In [20] it is shown that neglecting the Weyl 
focusing seems to be correct in our universe. 

The expansion 9 is related to the relative change of an 
infinitesimal surface area A of the beam's cross section 
by 



1 din ,4 

2 dX 



1 A 

2 A' 



(3) 



Since the angular diameter distance Da is proportional 
to y/A~, using Da instead of y/A and considering the affine 
parameter A as a function of the redshift z, it is obtained 
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T a pk a k?D A = 0, 

where we have used Einstein's equations to replace R a p 
in terms of T a p . 

A solution of the last expression is related to the an- 
gular diameter distance if it satisfies the initial condi- 
tions that make the wavefront satisfy Euclidean geometry 
when leaving the source: 
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where H is the present value of the Hubble parameter 
and c is the velocity of light. 

In order to take into account local inhomogeneities in 
the distribution of matter, Dyer and Roeder [10] con- 
sidered that the matter density p m should be inter- 
preted as the average density of the intergalactic regions 
where the observed light has traveled. Then to measure 
this effect they introduced a phenomenological parame- 
ter a = (1 — Pciumps// ) called the dumpiness or smooth- 
ness parameter, which quantifies the amount of matter in 



clumps relative to the amount of matter uniformly dis- 
tributed. Its maximum value a = 1 corresponds to the 
homogeneous FRW case, while for a totally clumped uni- 
verse a — 0; an intermediate value < a < 1 indicating 
a partially clumped universe. 

Assuming an energy-momentum tensor of the form 
T^ v = ap m u^u u + PAg^, where is the four- velocity 
of a comoving volume element, p m is the matter energy 
density and pa = A/87rG is the vacuum energy density 
associated to the cosmological constant, Eq. Q can be 
rewritten in the form 
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(6) 

We shall consider a flat (zero curvature) , homogeneous 
and isotropic geometry of Robertson- Walker (RW), in- 
cluding only dark matter (dust) and dark energy com- 
ponents. The corresponding Hubble parameter is given 
by 
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On the other hand, transforming the affine parameter 
to the observable redshift by using 
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and an appropriate change of variable, 

1 + n m z(3 + 3z + z 2 ) 
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h{D A ,z) = {l + z)D A , C(z) = 

(9) 

Eq. (|]) can be transformed into a hypergeometric equa- 
tion 

d 2 h fl 7 \ dh v(v + l) 



,36 



where the parameter v is related to the dumpiness pa- 
rameter, a, by 
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The range for v is < v < 2, where v — 0(a = 1) 
corresponds to a FRW fluid, while v = 2{a = 0) to a 
totally clumped universe. 

Since (1 — fi m ) < 1, the variable in Eq. ( |To| satis- 
fies C, > 1; thus to make it compatible with the defined 
range for the argument of the hypergeometric functions, 
we transform to the inverse argument, < 1, using 
properties of the hypergeometric functions. 

To connect the angular diameter distance Da with the 
observable luminosity distance e?^, we use the Ethering- 
ton relation [21], c?l = Da(1 + z) 2 . Incorporating also 
the initial conditions, the luminosity distance we used for 
the observational tests is given by [2"2"] . 
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FIG. 1: The plot corresponds to the ZKDR luminosity 
distance as a function of the redshift z, from Eq. (12); 



we have fixed fi m = 0.266 ± 0.029 from WMAP-7 years 
[39] and plot for different values of the smoothness 
parameter: a — (a completely clumped universe), 
a = 1 (homogeneous FRW) and for a partially clumped 
universe, a = 0.5. Clearly the effect of diminishing the 
smoothness parameter is to increase the luminosity 
distance. 



A few comments are in order: In Eq. ([lj, for the rate 
of expansion of the null congruence, the Ricci tensor, 
R a f3, makes null rays converge. Then it is reasonable to 
expect that diminishing R a p would result into reducing 
the convergence of the null bundle, i.e. there would be a 
dimming of light rays. This effect can be done by dimin- 
ishing the density of matter p m through the introduction 
of the smoothness parameter, p m t— > ap m , < a < 1. 
The plots in Fig. [T] for the ZKDR luminosity distance 
versus the redshift for different values of the smoothness 
parameter, show that this actually happens so: as a is 
smaller, grows faster. 

Have in mind that it is a naive reasoning, because nei- 
ther Eq.Q nor Eq.([3]) are linear equations. Notice that 
in this approach the luminosity distance, Eq. (12), does 



not totally capture the effect of inhomogeneities in the 
expansion rate, since in the derivation of cLl are assumed 
the same expressions for H(z) and X(z) than in the homo- 
geneous FRW model, Eqs. (|7|-Q. However, the effects 



of local inhomogeneities seem to be more complex to be 
described only through the smoothness parameter. 

In this regard, in Ref. 16 it is quantified the density 
probability distribution as a function of the beamwidth 
of supernovae and it is discussed a modified version of 
the DR approximation. It is shown that even for Gpc- 
lcnght beams of 500 kpc diameter, nonlinear corrections 
appear to be non-trivial, being not even clear whether 
underdense regions lead to dimming or brightening of 
sources (see Ref. (25). It was already noted by Linder 
in Ref. [53] that the smoothness parameter is a func- 
tion of the redshift. This dependence is also considerer 
in Ref. [25] , where a generalized Dyer-Roeder distance is 
introduced. However as was recognized in Ref. [26], this 
should be considered carefully because if the smoothness 
parameter changes with the redshift, the luminosity dis- 
tance expression will depend on several free parameters 
with difficult physical interpretation. A more natural de- 
pendence of the smoothness parameter with the redshift, 
coming from weak lensing, is presented in Ref. |26j . In 
Ref. [27] it was addressed the properties of cosmological 
distances in inhomogeneous quintessence cosmology. 

Although the DR approach has been criticized by sev- 
eral authors, the contrast of d L versus observational data 
so far has not lead to conclusive results, in the sense of 
totally discarding this model. Therefore it is important 
to elucidate its validity and scope in the description of 
the universe. So for the time being, we look for useful 
information from the simplest approach that considers a 
being constant. 



III. TESTING THE ZKDR DISTANCE VS. 
OBSERVATIONS 

A. Samples and statistical method 

We have used the updated observational data from 
Union2.1 supernovae data set reported in Ref. [13] and 
the distance modulus of long Gamma-ray Bursts re- 
ported in Ref. |14j to constrain the two free param- 
eters of the model, the matter density fi m , and the 
smoothness parameter v (related to a through Eq. ( 11 )). 
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In order to do that, we used a Markov Chain Monte 
Carlo (MCMC) code to maximize the likelihood func- 
tion C(6i) oc exp[— x 2 (#i)/2] where 0j is the set of model 
parameters and the expression for x 2 {@i) depends on 
the dataset used. The MCMC methods (completely de- 
scribed in Refs. [28H30] and references therein) are well- 
established techniques for constraining parameters from 
observational data. To test their convergence, here we 
follow the method developed and fully described in Ref. 
[31] . Additionally, we have tested the ZKDR luminosity 
distance with the observational Hubble parameter data 
given in Ref. [T7] . 



B. Supernovae Type la (SNe la) 

Hitherto we know one of the cosmological observations 
highly capable to measure directly the expansion rate of 
the universe are the Supernovae Type la. In this work, 
we use the updated compilation released by the Super- 
nova Cosmology Project (SCP): the Union2.1 compila- 
tion which consists of 580 SNe la [13]. The Union2.1 
compilation is the largest published and spectroscopically 
confirmed SNe la sample to date. Constraints from the 
SNe la data can be obtained by fitting the distance mod- 
uli fJ,(z). The distance modulus can be calculated as 



H(z j ) = 5]og w [D L (z j) 9 i )]+25 
= 5log w [d L (zj,6i)] + no, 



(13) 



where /io = 42.38 — 51og 10 /i, h is the Hubble constant 
Ho in units of 100 km s _1 Mpc _1 , and diXzj,9i) is the 
ZKDR luminosity distance given by Eq. ( fl2| ); 9i denotes 
the vector model parameters (Q m ,z/). The x 2 function 
for the SNe la data is given by 
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where the a 2 corresponds to the error on distance mod- 
ulus for each supernova. The parameter fj,Q m Eq. ( 13 ) 



is a nuisance parameter since it encodes the Hubble pa- 
rameter and the absolute magnitude M, so it is more 
convenient to marginalize over it. 

Here, we maximize the likelihood by minimizing \ 2 
with respect to /«o, as suggested in j3U[33]. Then one 
can rewrite the \ 2 function as 



XsnW) = ci - 2c 2 Mo + C3/U0, 
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The minimization over no gives [1q — C2/C3. So the x 2 
function takes the form 



XsAOi) = ci - ~f. 

C3 



(19) 



Since Xsn = Xsw(Mo = 0, ftj) (up to a constant), we 
can instead minimize Xsn which is independent of fig. 



C. Long gamma-ray bursts (LGRBs) 

Previous analysis have shown that the ZKDR distance 
at short distances (redshifts in the range 0.1 < z < 1.7) 
does not make a difference respect to the ACDM model; 
it is then appealing to test the model for larger distances, 
like the ones of the GRBs, that extend the redshift range 
as far as 8.1. However, the problem is that GRBs appear 
not to be standard candles and to extract cosmological 
information it is necessary calibrate them for each cos- 
mological model tested. 

There have been several efforts to calibrate the corre- 
lations between the luminosity and spectral properties of 
GRBs in a cosmology-independent way and some propos- 
als to use SNe la measurements to calibrate them exter- 
nally are given in Refs. [54TI36] . 

More recently, in Ref. 14 , it has been estimated the 
distance modulus to long gamma-ray bursts (LGRBs) us- 
ing the Type I Fundamental Plane, a correlation between 
the spectral peak energy E p , the peak luminosity L p , 
and the luminosity time Tl = Ei so /L p , where Ei so is the 
isotropic energy. Basically the calibration was done in 
this way: first, the Type I Fundamental Plane of LGRBs 
was calibrated using 8 LGRBs with redshift z < 1.4 and 
SNe la (Union2) in the same redshift range by a local 
regression method, to avoid any assumption on a cosmo- 
logical model; then this calibrated Type I Fundamental 
Plane was used to measure the distance modulus to 9 
high-rcdshift LGRBs (see Ref. [H] for calibration's de- 
tails). We used the 9 calibrated LGRBs reported in Ta- 
ble 2 of Ref. [13] to derive constraints on tt m and v from 
the ZKDR luminosity distance. The \ 2 function for the 
GRBs data is defined by 
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where n{zj) = 51og 10 [di(2^, 9i)] + 25. Notice that we 
have used the standard expression of x 2 given through 
the observed distance moduli just to be consistent with 
the way in which the calibration was done. We have fixed 
H Q as 73.8 ± 2.4 from Ref. [37]. 
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D. Hubble parameter observations 

We use the compilation of Hubble parameter measure- 
ments estimated with the differential evolution of pas- 
sively evolving early-type galaxies as cosmic chronome- 
ters, in the redshift range 0.09 < z < 1.75 recently up- 
dated in Ref. [17] but first reported in Ref. [38] , 

This approach consists in the measurement of the dif- 
ferential age evolution of these chronometers as a func- 
tion of redshift, in this way a direct estimate of the 
Hubble parameter H(z) = -1/(1 + z)dz/dt ~ -1/(1 + 
z)Az/At is obtained, providing a reliable differential 
quantity Az/At, with many advantages in minimizing 
common issues and systematic effects. Compared with 
other techniques, it provides a direct measurement of the 
Hubble parameter, and not of its integral, in contrast to 
SNe la or angular/angle-averaged BAO. 

Observed values of H{z) can be used to estimate the 
free parameters of the model and also the best-value for 
H by minimizing the quantity 



x%(H ,e i ) = 



i) — H b s (Zj)] 
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where are the measurement variances. Hq has been 
fixed from Ref. [57]. ifn = 73.8 ± 2.4. The vector of 
model parameters, 0i, in our case will be (9; = (f2 m , v). 

Noteworthy that constraining the free parameters of 
this model via H(z) is controversial because this is not 
a direct comparison with the ZKDR luminosity distance, 
but we have calculated H(z) from the expression 
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connecting in this way H(z) with the parameter v (or a). 
We note that we are mixing the homogeneous FRW H(z) 
with the inhomogeneous proposal, however a joint analy- 
sis with other probes makes sense and has been addressed 
in previous works [18 . 



IV. RESULTS 



The results of our analysis are shown in Tables [T| |VIII| 
We made the analysis with the parameter i/, but in sev- 
eral tables are included the results for both parameters, 
v used by Kantowski, and a related with v by Eq.( 111. 

In Tables [T| and [H] are displayed the best fits for S2 m , 
v (or a) from the SNe la and GRBs probes reported 
in Ref. [13] and Ref. [H], respectively. Table [I] cor- 
responds to the analysis assuming a Gaussian prior on 
fl m = 0.266 ± 0.029 from WMAP-7 years [33] ■ Fig 
shows the corresponding contour levels in the (Q m , 
plane. 

Table [H] is the same as Table [I] but without assuming 
any prior. The inclusion of a prior does not have a dra- 
matic effect, but there is an improvement in the GRBs 



Sample 
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SNe la 
GRBs 
Joint 


n 98 ,+ - 019 

n ocq+0.028 
u.zj»_ 02g 

n 9SzL+ 021 


n w+ ' 417 
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0.856+° .f 7 6 6 
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0.975 
0.877 
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TABLE I: Summary of the best estimates of model 
parameters (Q m , v), obtained from the ZKDR 
luminosity distance using a prior on Sl m . The respective 
samples are SNe la reported in Ref. [13] and GRBs 
reported in Ref. [14] . The errors are at 68.3% 
confidence level. Joint stands for the joint analysis SNe 
la + GRBs. The corresponding confidence regions are 
shown in Figure [2] 
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0.975 
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TABLE II: Summary of the best estimates of model 
parameter (fi m , v). Same samples as in Table |l| but, m 
this case, no prior in f2 m is assumed. The errors are at 
68.3% confidence level. 



analysis, that goes from a Xrcd = 0.84 without prior to 
Xred — 0.877 when the prior is assumed. The values 
of the smoothness parameter show a degree of dumpi- 
ness in the cosmic fluid, however it does not have an 
impact on reducing the proportion of dark energy, that 
remains over the 70% of the total density; more precisely 
Q\ = 1 — Q m = 0.7161q q2i from the joint analysis of 
SNe la and GRBs. These results make us think that the 
DR approach does not describe the effects of dumpiness 
in a reliable way. 

Tables [III] and |TV] correspond to the 2a confidence level 
constraints of the cosmological parameters, Q rn and a 
with and without the prior on fl m , respectively. 

Table [TV] showing the 2-er confidence level for fi TO and 
a without assuming a prior on Sl m , can be directly com- 
pared with the results given in Ref. [TS] (that we include 
for completeness in table Ivb. In Ref. [15J it is confronted 
the ZKDR luminosity distance versus SNe la and GRBs 
probes with the samples data 557 SNe la Union2 [3D] and 
the 59 Hymnium GRBs 41j. Using a x 2 minimization, 
their results, shown in Table [V] did not constrain the 
smoothness parameter. The best fits obtained by com- 
bining the SNe la and GRBs probes are f2 m = 0.27 and 
a = 1, with a x 2 = 568.36 or X 2 cd = 0.927. 

Particularly, for GRBs, using a better sample and the 
MCMC method, we obtained a reliable constraint with 
a x 2 cd = 0.877 from GRBs only, and X 2 cd = 0.973 from 
the joint analysis without assuming the prior on £l m and 
X 2 c d = 0.975 from the joint analysis assuming the prior 
on fl m . 

In a previous work done by Busti et al. |18) . a chi- 
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FIG. 2: Confidence regions in the (fi m , v) plane for the 
model with a ZKDR luminosity distance using a prior 
on tt m . The contours correspond to \a-2a confidence 
regions using: LGRBs, largest region on the back; SNe 
la, smallest region on the front; the combination of the 
two observational data, the region between the LGRBs 
region and SNe la region. 
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f2,„(2 
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Xred 


SNe la 


0.250 < o m 


< 0.325 


0.461 < a < 0.995 


0.975 


GRBs 


0.204 < fi m 


< 0.315 


0.183 < a < 0.952 


0.877 


Joint 


0.246 < n m 


< 0.327 


0.344 < a < 0.965 


0.975 



TABLE III: Limits to fl m and a with prior on f2 m . The 
errors are at 95% confidence level. Joint means the 
combined analysis with SNe la + GRBs. 



Sample gmgtr) a{2a) Xred 

SNe la 0.252 < Q. m < 0.373 0.247 <a< 0.991 0.975 
GRBs 0.142 <fl m < 0.447 0.072 < a < 0.988 0.840 
Joint 0.248 < fl m < 0.385 0.150 < a < 0.940 0.973 



TABLE IV: Limits to Q m and a without prior on fl 
The errors are at 95% confidence level. 



Sample Q m (2a) a(2a) v(2a) x?od 

SNIa 0.24 < n m < 0.35 0.25 < a < 1 0<v< 1.68 0.984 
GRBs 0.11 < Q m < 1 unconstrained unconstrained 0.411 
Joint 0.24 < Q, m < 0.33 0.52 <a<l 0<^< 1.27 0.927 



TABLE V: Obtained constraints for (f2 m ,a) in Rcf. 
[15] ; SNe la data from Ref. [10] and GRBs data taken 
from Ref. |41j . The smoothness parameter remains 
unconstrained from GRBs. 



On (2a) a(2a) xled 

No Prior 0.208 < Q. m < 0.376 0.510 < a < 0.995 1.032 
With Prior 0.223 < Q. m < 0.315 0.626 < a < 0.996 1.025 



TABLE VI: Limits to O m and a from H(z). The 

second row is obtained assuming a prior of 
n m = 0.266 ± 0.029 from WMAP-7 years [39]. The 
errors are at 95% confidence level. 



Sample Q, m ct 

182 SNe la +Hubble 0.27 < fi m < 0.37 0.66 < a < 1 
12 Hubble 0.27 < Q m < 0.3 0.5 < a < 1 



TABLE VII: Summary of the best estimates of model 
parameters for the Hubble probe comparing the 
previous results by Busti (2011) [18] from samples 
reported in Ref. [12] and Ref. [13] and in the second 
row the ones presented in here only from the Hubble 
parameter given in the sample by Moresco (2012) [T7] 



square analysis was performed with Hubble parameter 
measurements reported in Ref. [42] and in Ref. [43] . Our 
analysis provides better constraints, that are summarized 



in Table [VT] while in Table VII we make the comparison 
between the results obtained in Rcf. [18] and ours. 

Considering that the samples of SNe la and GRBs cor- 
respond to disjoint intervals of redshift, the different val- 
ues of a for each probe can be interpreted as a depen- 
dence of the dumpiness parameter a on the redshift (see 
Table |VHl|. 



Sample Q m a redshift range Xred 

SNe la 0.2851^1 0.856tg j?| 0.015 < z < 1.414 0.975 

GRBs 0.259±g 0-587±gJ^ 1.547 < z < 3.57 0.877 

Hubble 0.268tg;^3 0.895^ ^| 0.09 < z < 1.75 1.025 

Joint 0.275to;ois 0.821^129 0.0 15 < z < 3.57 0.974 



TABLE VIII: Summary of the best estimates of model 
parameters and the corresponding redshift range using 
in all the cases a prior on 57 m from Ref. [39]. The 
smoothness parameter a shows a dependence on the 
redshift range. Joint: SNe la + Hubble + GRB. 
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V. CONCLUSIONS 

We probed the ZKDR luminosity distance proposed 
by Dyer and Roeder (DR), that includes the effect of lo- 
cal inhomogeneities through the smoothness parameter 
a, with data from SNe la [T3] and GRBs [HJ. Addition- 
aly we tested with the direct Hubble measurements 
reported in Ref. [T?] . 

Better data and a refined statistics for GRBs allow 
us to probe the ZKDR luminosty distance in a range 
of higher rcdshifts with confident results. Indeed, our 
analysis improves, in great deal, previous works that did 
not succeed in constraining the smoothness parameter. 

We can summarize the presented results as follows: 

1) We have probed in a reliable way two distinct red- 
shift ranges, in fact two disjoint intervals of z: the SNe 
la range, 0.015 < z < 1.414 and 1.547 < z < 3.57 for 
GRBs. 

The obtained results, two different values of a in each 
interval, strongly suggest a dependence of the dumpiness 
parameter on the redshift, see Table |VIH| for the values 
of a and the corresponding redshift range. Since the red- 
shift intervals are disjoint, the different values of a show 
a kind of evolution in the smoothness parameter with 
respect to the redshift. In principle, assuming that this 
approach describes the dumpiness effects with accuracy, 
one would expect that nearby regions should correspond 
to more dumpiness, then to smaller values of a. In other 
words, the structure formation process leads to a more lo- 
cally inhomogeneous universe and thus, the smoothness 
parameter should evolve from homogeneity (a = 1) to 
total dumpiness (a = 0). However, the test turns out 
on the contrary: the smoothness parameter evolves to- 
ward homogeneity, being a — 0.587T^202 f° r the range 
1.547 < z < 3.57 and a greater value, a = 0.856^^^ for 
the SNe la range, 0.015 < z < 1.414. These results show 
that the Dyer-Roeder equation for should be revised, 
as several authors have suggested [It)] . 



2) Our results indicate the existence of dumpiness, be- 
cause the values of a differ from the homogeneous FRW 
model (a = 0.856, a = 0.587 and a = 0.895 from SNe la, 
GRBs and Hubble, respectively), however this does not 
have a clear impact in the amount of dark energy needed 
to fit the observations. According to SNe la, Hubble and 
GRBs probes, the amounts of dark energy that fit to ob- 
servations are: Q A = 0.715^;^, fl A = 0.732 ± 0.023 
and Ha — 0.741 ± 0.028, respectively (with prior on 
51 m ); in other words, the average matter density in- 
creases in the range of small z (Q m = 0.285^q'oi8i fr° m 
SNe la compared with the value in the range of GRBs, 

i£ m — o.028/ ' 

In conclusion, our results indicate that DR approach 
models in an incorrect way the effect of local inhomo- 
geneities in the propagation of light. In Section II we 
have noticed that the model is formulated with at least 
two questionable assumptions: that the Hubble parame- 
ter and the affine parameter as a function of the redshift 
have the same expressions than in the homogeneous FRW 
model. 

So far, the evidence indicates that a modification of 
the Dyer-Roeder equation is necessary to describe in a 
reliable way the backreaction effects due to local inho- 
mogeneities. It is worth then to test with observational 
data some of the improved models that try to mend the 
drawbacks of the DR model; for instance in Ref. [Mj it is 
considered a correction in the analogous of the smooth- 
ness parameter coming from density fluctuations along 
the line of sight. 
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